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Abstract Avian infectious bronchitis virus (AIBV) is 
classified as a member of the genus coronavirus in the family 
coronaviridae. The enveloped virus has a positive-sense, sin- 
gle-stranded RNA genome of approximately 28 kilo-bases, 
which has a 5’ cap structure and 3’ polyadenylation tract. 
The complete genome sequence of infectious bronchitis virus 
(IBV), Beijing isolate, was determined by cloning sequencing 
and primer walking. The whole genome is 27733 nucleotides 
in length, has ten open reading frames: 5'-orfla-orflab-s-3a- 
3b-e-m- 6a-6b-n-3'. Alignments of the genome sequence of 
IBV Beijing isolate with those of two AIBV strains and one 
SARS coronavirus were performed respectively. The genome 
sequence of IBV Beijing isolate compared with that of the 
IBV strain LX4 (uncompleted, 19440 bp in size) was 91.2% 
similarity. However, the full-length genome sequence of IBV 
Beijing isolate was 85.2% identity to that of IBV Strain 
Beaudette, and was only 50.8% homology to that of SARS 
coronavirus. The results showed that the genome of IBV has 
remarkable variation. And IBV Beijing isolate is not closely 
related to SARS coronavirus. Phylogenetic analyses based on 
the whole genome sequence, S protein, M protein and N pro- 
tein, also showed that AIBV Beijing isolate is lone virus in 
group [II and is distant from SARS coronavirus. In conclu- 
sion, this study will contribute to the studies of diagnosis and 
diseases control on IBV in China. 
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Avian infectious bronchitis virus (AIBV) is a highly 
infectious, contagious pathogen of chicken and causes 
disease of the respiration and urogenital tracts, resulting in 
significant economic losses to poultry industries through- 
out the world!"!, In China, IBV leads to more than RMB 
one billion yuan losses every year. IBV, as a member of 
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the family coronaviridae in the new order Nidovirales, has 
a single-stranded positive-sense RNA genome, which is 
27.6 kb in length’!. The genome of IBV has more than 10 
open reading frames, which encode two polyproteins and 
three major structural proteins and other small proteins. 
Three major structural proteins are the spike glycoprotein 
(S), the membrane glycoprotein (M), and the phosphory- 
lated nucleocapsid protein (N). The S protein is cleaved to 
release S1 and S2 proteins after translation. The S1 protein 
induces protective immunity’ *!, And the most antigenic 
of the virus neutralization antibody-inducing epitopes is 
formed by a few amino acids that occur in the first and 
third quarters of the S1 subunit, especially in the hyper- 
variable regions!*”!, 

IBV is the lone virus in the antigenic group III of 
coronaviridae'*!, however, more than 20 serotypes within 
IBV have been recognized worldwide on the basis of vi- 
rus-neutralization test!!. Therefore, it is very difficult to 
prevent from the infections of virus. The whole genome of 
IBV was determined in 1987”!. And the genome se- 
quences of IBV are very important and useful for disease 
diagnosis and vaccination. However, remarkable variation 
was observed in various IBV isolates: the genome of IBV 
isolate in this region may be very different with that in the 
other regions. IBV Beijing isolate, which caused great 
economic losses, was obtained from a chicken house in 
Yukou, Beijing!'”!, 

In additional, an outbreak of atypical pneumonia, 
named severe acute respiratory syndrome (SARS) by the 
WHO, was identified in Guangdong Province, China, and 
spread to several countries. This disease is great threaten- 
ing to the health of people worldwide. Many reports 
showed that a novel coronavirus was the pathogen of the 
SARS"''’!_ And scientists considered that SARS coro- 
navirus could be from some animals. After the complete 
genome of IBV Beijing isolate was sequenced, we deter- 
mined the relations between IBV Beijing isolate and 
SARS coronavirus. 


1 Materials and methods 


(i) Virus. IBV Beijing isolate''®! was propagated 
in 10-d-old specific-pathogen-free (SPF) embryonated 
chicken eggs. The allantoic fluid was collected 24—36 h 
after inoculation, clarified by centrifugation at 5000 r/min 
for 30 min, 4°C. The pellet was discarded and the super- 
natant was concentrated by more than 2 h high-speed cen- 
trifugation (30000 1/m) at 4°C. The supernatant was dis- 
carded and the pellet was resuspended in PBS (pH 7.2). 
The virus suspension was ultra-centrifuged through a 
30%/40%/50% saccharose step gradient at 32000 r/min 
for 3 h at 4°C. Mixture between 40% and 50% saccharose 
was transferred and PBS (pH 7.2) was added. Samples 
were mixed and concentrated by centrifugation at 30000 
r/min for 2 h at 4°C. The supernatant was discarded and 
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the virus pellet was resuspended in PBS. The virus sus- 
pension was divided into 250 uL each tube and was stored 
at -80°C_ until use. 

(ii) Viral RNA extraction and RT-PCR. Viral ge- 
nomic RNA was extracted with Trizol LS reagent (Invi- 
trogen, USA). Each 250 pL virus suspension was added 
with 750 uwL Trizol LS reagent and mixed. The homoge- 
nized samples was incubated for 5 min at 15°C—30°C 
and then added with 200 wL chloroform. The capped sam- 
ple tubes were shaken vigorously by hand for 15 s and 
incubated at 15°C—30°C for 2 min. The samples were 
centrifuged at 1100 r/min for 15 min at 4°C. The upper 
aqueous phase was transferred to a clean tube and added 
with the same volume of isopropyl alcohol. The samples 
were incubated at 15°C—30°C for 10 min after being 
mixed. The mixture was centrifuged at 1100 r/min for 10 
min at 4°C. The supernatant was discarded and RNA pellet 
was washed twice with 75% ethanol. At the end of the 
procedure, the RNA pellet was air-dried (about 15 min) 
and dissolved by 10 uL 0.1% DEPC-treated water. 

The random hexameric primer (Promega, USA) or a 
specific reverse primer was used in cDNA first strand 
synthesis. After cDNA first-strand was synthesized, each 
specific primer pair was used in the reaction of PCR. A 
20-uL reaction volume was used for the first strand cDNA 
synthesis. Total RNA (10 uwL) and random hexameric 


primer (0.5 ug/uL, Promega, USA) or a specific reverse 
prime 1 pL were added to a nuclease-free microcentrifuge 
tube. Mixture was heated to 70°C for 10 min and was 
quickly chilled on ice for 1—2 min. After brief centrifuga- 
tion, the following components: 4 uL 5 *X first-strand 
buffer, 2 uL 0.1 mmol/L DTT, 1 pL RNaseOut (40 U/pL) 
and 1 uL M-MLV reverse transcriptase (200 U/uL) were 
added. After being mixed gently, the mixture was incu- 
bated at 37°C for lh. Samples were inactivated at 70°C 
for 15 min and subsequently treated with 2 U of Ribonu- 
clease H at 37°C for 30 min. The cDNA was used as tem- 
plate for amplification in PCR. A 30 uL reaction volume 
was used in PCR. The components were 1.5 uL cDNA, 3 
uL 10 X buffer, 2 uL 2.5 mmol/L dNTP each, | uL of each 
primer (25 pmol/L), 0.5 uwL Taq DNA polymerase (5 U/uL) 
and 21 uL ddH,0. PCR was carried out with the following 
thermocycling parameters: denaturation at 94°C for 5 min, 
35 or 40 cycles of denaturation at 94°C for 40 s, anneal- 
ing at 49°C—53°C for | min, and extension at 72°C for 1 
—2 min or 68°C for 4—5 min, a final extension at 72°C 
for 10 min and store at 4°C. Annealing and extension con- 
ditions were dependent on primers and the length of 
amplification fragments. Some primers used in PCR are 
listed in Table 1. 

Rapid amplification of cDNA ends (RACE) was 
performed to capture the 3’ and 5’ ends of the viral ge- 


Table 1 Primers for PCR 
Primer Position/nt Size/nt Forward primer Reverse primer 
gpst 1—355 355 ACT TAA GAT AGA TAT TAA TAT A GGT GAC CCC GTA/G TTT C 
gpl 10421928 887 CAT GGA TGC TT/CT GGG CTC AAA T/C TCC/T TCA AGA/G TTT/G CGC ACT TTT CC 
gp3 28176014 3198 =G/ACG GTAAAG AC/TT TGG ATT GC/TA T GCA/G TGT TTA/G TAC AAA TGA AGC 
gp3-2 27933540 748 AGAATGTC/TG CAC TTG TTG ATAA TT/CT TAG CC/TT TCT/C TCT CAC GTT 
gp3-3 40506014 1965 TGAG/ACC AAG AA/GC AC/TA/G TCG ATTA = AA/GG CG/AT GTT TG/AT ACA AAT GAA G 
gp4 6708— 10454 3747 TTC TTAAGGAT/AG CAC A/TA/GA AAT G TAA CA/CA TG/AA CA/TG CCA TTA GA/TA C 
gp4-2 69609601 2642 ATT ATAAGG CA/GG GCA CAC TTC G ATT G/ATC A/CC/G CA/TG CCC ACT TAT TA 
gp5 10498 — 14242 3745 = TTIGTIACTT TGGA/GTT C/GTA TTGG CTA TAC GCA ACA AAT TT/AG GCA T 
gp5-2 12341 12956 616 TCT GTT CAA TCA GTT GCT GGA/T CCA ACA TAA CCT/C TTT TCA ACCA 
gpSa 1459417221 2627 AAAAGTTTTATT CTT ACA/T TG/AT G CTT CTA CAT CAA AAC CNA C 
gpSa-2 14896— 16710 1815 GTAGAT GAC/A GTG GAT AAG ACA CTG G/ATT CAT AGC ATT GTA AGG T 
gp6 1761218097 486 TTAATT CTCATA CTC AG/AG A/CTTA ACA TTG/A CAA TTA/C CAA AAC ATA C 
gps 1844221319 2878 = TAT/C AAT GCA GCT GTT/C ACA TAG/A AAT TAT AAT AAC CAC TCT 
gp9 2170623182 1477. =TTTAG/AA TAAI/C G/TTG TGT TGA G/TTA T GTA TGG TTA A/GG/AA CAT GTC T/GT 
gp9-3 2219322681 489  TTC/T AAC/T TTA/G ACT GTT ACA GAT G CCA TAG AA/TG CTA CTA G/AAG AA 
gp3t-2 2518425901 718 GAT ACT GGC GAG CTA GAA AG ACG CGG GAC AAG GCT CTG 
gp3t-3 2620927445 1237. GGT TTAAGC CAG GCA AAG GT CTG GCG ATA GAC ATG TAG GT 


nome. A 20-uL reaction volume was used for the first 
strand cDNA synthesis in 3’ RACE. The procedure was 
done as follows: AP (5’-GGCCACGCGTCGACTAGTA- 
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CTTTTTTTTTITTTTTTTT-3') 1 pL, mixture contained 
10 pL total RNA, were heated to 70°C for 10 min and 
was quickly chilled on ice for 1—2 min. After brief cen- 
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trifugation, 2 uL 10 buffer, 2 wL 0.1 mol/L DTT, 2 ul 
25 mmol/L MgCL,, | wL 10mmol/L dNTP Mix, 1 ul 
RNaseOut (40 U/uL) and 1 pL M-MLV reverse transcrip- 
tase (200 U/L) were added. Mixture was mixed gently 
and incubated at 37°C for 1 h. Samples were inactivated 
at 70°C for 15 min and subsequently treated with 1 pL 
Ribonuclease H (2 U/uL) at 37°C for 30 min. 1.5 uL 
cDNA was used as template in a 25-uL PCR reaction 
volume. 0.5 pL AP, 0.5 pL GSP (5'-ATTACCTACATGT- 
CTATCGCCA-3'), 2.5 uwL 10 X buffer, 25 mmol/L MgCL 
1.5 pL, 0.5 wl 10 mmol/L dNTP Mix, Taq DNA poly- 
merase (5 U/L) 0.5 wL, 17.5 pL ddH,O were mixed gen- 
tly and heated to 94°C for 3 min to denauture. 35 cycles 
were carried out for 30 s at 94°C, 1 min at 62°C, 1 min at 
72°C, with a final extension at 72°C for 10 min, then 
stored at 4°C. 

A 25-uL reaction volume was used in 5'RACE for 
the first strand cDNA synthesis. Mixture consisted of 15.2 
uL total RNA, 0.3nL GSP1 (5'-ACTTAAGATAGATA- 
TTAATATA-3', 25 umol/L) was heated to 70°C for 10 
min and quickly chilled on ice for 1—2 min. After brief 
centrifugation , 2.5 uL 10 X buffer, 2.5 uL 0.1 mol/L DTT, 
2.5 uwL 25 mmol/L MgCL», | wL 10 mmol/L dNTP Mix, 1 
uL SuperScript I] RNase H-RT (200 U/L) were added. 
Mixture was mixed gently and incubated at 42°C for 1 h. 
Samples were inactivated at 70°C for 15 min and subse- 
quently treated with 1 uL RNase Mix (2 U/uL) at 37°C 
for 30 min. cDNA was purified and tailed before PCR. 
The two steps were done according to the instruction 
manual. A 30 uL reaction volume contained 3 pL cDNA, 
3 uL 10Xbuffer, 2 wl 2.5 mmol/L dNTP each, GSP2 
(10 umol/L, 5'-CGCTACTTCCACACACCC-3') 0.5 uL, 
AAP(10 umol/L, 5'-GGCCACGCGTCGACTAGTACG- 
GGIIGGGIIGGGIIG-3") 0.5 wL, 0.4 uL Taq DNA poly- 
merase (5 U/uL) and 21.3 pL ddH,O was used for ampli- 
fication in PCR. PCR was carried out with the following 
thermocycling parameters: denaturation at 94°C for 3 min, 
35 cycles of denaturation at 94°C for 30 s, annealing at 
57°C for 45 s, and extension at 72°C for 1 min, a final 
extension at 72°C for 7 min and stored at 4°C. 

(iii) Sequencing and analysis. PCR products were 
cut from 1% agarose gels and purified using the GENE- 
CLEAN III Kit (Qbiogene). The purified PCR products 
were either sequenced directly or cloned into the PMD 
18-T vector (TaKaRa Biotech. Co., Itd., China), and auto- 
mated sequencing with BigDye Terminator kit (Perkin 
Elmer) run on an ABI 377 DNA Sequencer. Primer walk- 
ing is the most common approach for directed sequenc- 
ing!'*!, and was used mostly in this study. 

Initially, eleven primer pairs to various regions of the 
genome were designed according to the full-length ge- 
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nome of IBV strain Beaudette (Accession: NC_001451.1) 
and the PCR products were 1.5 kb in size. Many products 
were obtained. However, none was the target fragments. 
The PCR products were purified and cloned into the PMD 
18-T vector before sequencing. Each clone was sequenced 
using M13F or M13R. 160 sequences were obtained and 
used to search the related sequences in the GenBank using 
the BLAST software  (http://www.ncbi.nlm.nih.gov/ 
BLAST). The viral sequences were kept and compared to 
the complete genome of IBV strain Beaudette (Accession: 
NC_001451.1). Positions of the sequences were deter- 
mined to have high similarity more than 80%. Twelve 
contigs and 13 gaps were got after these sequences assem- 
bled. According to these sequences, we designed addition- 
ally 12 primer pairs for gaps filling. The largest gap was 
less than 4 kb. Each primer pair was used for PCR and had 
a target fragment. These fragments were purified and both 
strands were sequenced. Sequences of gaps were deter- 
mined by primer walking. Other 77 primers were designed 
for PCR and sequencing, but not all were used for PCR. 
Finally, the complete sequence was assembled, by using 
Phrap (http://www.phrap.org). 


2 Results 


The complete genome sequence of IBV Beijing iso- 
late was determined using PCR and cloning before se- 
quencing as well as primer walking approaches. After 
clones (recombinant plasmids) were sequenced, the as- 
sembly of 140 sequences were totally 11 kb in length. 
These sequences comprised 12 contigs, covering these 
regions more than 4-fold. 101 primers were used for direct 
sequencing, among which 36 primer pairs were used for 
PCR and 3 primers for RACE. 136 sequences were ob- 
tained by primer walking. The sequences assembled were 
about 17 kb in size. Totally, 276 sequences resulted in the 
complete 27733 bp sequence of IBV Beijing isolate ge- 
nome. The average useful length of the component se- 
quences was 350 bases. The IBV Beijing isolate genomic 
sequence has been deposited into GenBank (Accession: 
AY319651). 

(i) Genome organization. The 3’ and 5’ ends of 
the viral genome were captured after rapid amplification 
of cDNA ends (RACE) was performed. A fragment about 
300 bp in size from 3’RACE and a band about 350-bp- 
length from 5'RACE were purified and sequenced directly 
or cloned into the PMD 18-T vector before sequencing. 
The sequences of RACE were almost 300 bp in length, 
with a 17polyA tract in 3’RACE, and 13polyG tail in 5’ 
RACE (Fig. 1). 

The IBV genome has five major genes, with a typical 
gene order: rep, spike, envelope, membrane, nucleocapsid. 
The rep gene, which comprises approximately two-thirds 
of the genome, is encoded two polyproteins (encoded by 
orfla and orflb). The four genes downstream of rep en- 
coded the structural proteins S, E, M and N. There are also 
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ATTACCTACATGTCTATCGCCAGGGAAATGTCTAATCTGTCTACTTAGTA = 50 

51 GCCTGGAAACGAACGGTAGACCCTTAGATTTTAATTTAGTTTAATTTTTA 100 
101 GTTTAGTTTAAGTTAGTTTAGAGTAGGTATAAAGATGCCAGTGCCGGGGC 150 
15] CACGCGGAGTACGATCGAGGGTACAGCACTAGGACGCCCACTAGGGGAAG 200 
201 AGCTAAATTTTAGTTTAAGTTAAGTTTAATTGGCTAAGTATAGTTAAAAT —- 250 
251 TTATAGGCTAGTATAGAGTTAGAGCAAAAAAAAAAAAAAAAA 


(a) 


GCACTCGAACCGACCACGGCGTCCGCAGGACGCTGGCACTAGGTAGCCCA 50 
51 AGCAACGGATGTAATAAAGTGAGACACAACAAAACAAGCAGTGATACTAG — 190 
101 CAACAATAAAAATTTAATAACTAAAACCTGACAGGTGGCCAGGTGCCATC = 150, 


151 CAGGGCACTAAAGTGCACTGCAATGGAACACCTATAGGGACCAGCTGTAG 200 
201 GTATTTAAGTCCGTTTTGTTAAGTTGGAAATCTAGCGCAAGGCTAGTGCA 250 
251 ATAGATATATATTAATATCTATCTTAAGTCCCCCCCCCCCCC 


(b) 
Fig. 1. Sequences obtained by RACE. (a) Sequence of 3'RACE; (b) 
sequence of 5'RACE. 


= 


genes encoding a number of nonstructural proteins that are 
located between S and E, and between M and N. The ge- 
nome organization of IBV Beijing isolate was determined 
and ten open reading frames (ORFs) were predicted, in- 
cluding genes encoded unknown proteins: 3a, 3b, 6a and 
6b (Table 2, Fig. 2). Other ORFs encoded unknown small 
proteins were presumed (data not shown). 

(ii) Comparison of sequences. To evaluate how 
conservative the IBV genome was and how homologous 
the sequences were between IBV Beijing isolate and 
SARS coronavirus, the sequences of two strains of IBV 
and one SARS coronavirus were downloaded from Gen- 
Bank. The complete genome of IBV Beijing isolate com- 
pared with that of IBV strain Beaudette (Accession: 
NC_001451.1) resulted in a similarity of 85.2%. However, 
there was 91.2% identity between IBV Beijing isolate and 
IBV LX4 (Accession: AY223860), whose genome had not 
been sequenced completely and was only 19940 bp in 


ORF la 


size. And only 50.8% homology was shown between the 
full-length genome of IBV Beijing isolate and that of 
SARS coronavirus (NC_004718.1, TOR2). 


Table 2 Ten ORF of IBV Beijing isolate 


ORF Position/nt Peptide length/a.a. 
530—12391 
RNA polymerase 12391—20454 6641 
Orfla 530—12421 3963 
S protein 20405—23914 1169 
3a protein 23914—24087 57 
3b protein 24087—24275 62 
E protein 24259—24588 109 
M protein 24557—25237 226 
6a protein 25601—25798 65 
6b protein 25795—26040 81 
N protein 25983—27212 409 


Sequences encoding five proteins of IBV Beijing 
isolate were also compared with those of IBV strain 
Beaudette and SARS coronavirus, at the levels of nucleo- 
tide and amino acid, respectively (Tables 3 and 4). S pro- 
tein and E protein had low identity at the levels of nucleo- 
tide and amino acid between two strains of IBV. While 
RNA polymerase was of relatively high conversation at 
the amino acid level between two strains. However, four 
structural proteins of IBV Beijing isolate had fairly low 
homologies (<50%) to those of SARS coronavirus. On the 
contrary, the homologies of four structural proteins at the 
amino acid level were lower than those at the nucleotide 
level, respectively, between IBV Beijing isolate and SARS 
coronavirus. 
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Genome organization of IBV Beijing isolate. Overall organization of the IBV Beijing isolate is 27733 nt in length. Predicted ORFs la and lab, 


encoding the polyproteins, and those in turn encoding the S, E, M and N structural proteins are indicated. ORFs encoding unknown proteins 3a, 3b, 6a, 
6b are also shown. The vertical position of the boxes indicates the phase of the reading frame, and the known coding regions (grey boxes) and reading 


frames encoding potential proteins (oblique boxes) are indicated. 
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Table 3__nt identity of five proteins (%) and three structural proteins of IBV Beijing isolate with 

Gene AIBV” SARS-CoV” those of other animal coronavirus and SARS coronavirus 
orfla 83.5 48.9 for which sequences information was available in the 
orflb 88.7 59.2 GenBank. The topologies of the resulting phylogenetic 
spike 81.5 49.8 trees are remarkably similar (Fig. 3). In all cases, SARS- 
envelope 85.5 48.0 CoV and IBV Beijing isolate are not closely related to 
membrane 87.5 49.0 other known coronavirus. IBV Beijing isolate is a lone 
guclee capsid Be pau virus in group III"! and SARS was segregated into a 


a) AIBV strain Beaudette (Accession: NC_001451.1); b) SARS 


coronavirus (Accession: NC_004718.1/TOR2). fourth branch. 


Table 4 _a.a identity of five proteins (%) 3 Discussion 


Protein AIBV” SARS-CoV” There is very remarkable variation in the genome 
Orfla 91.1 31.5 sequence of IBV. The full-length genome of IBV Beijing 
Orflb 96.8 67.2 isolate, which was 27733 bp in length, was 125 bp longer 
Spike 86.5 32.4 than that of IBV strain Beaudette (Accession: NC_ 
Envelope 86.2 29.1 001451.1), sharing only 85.2% identity of the nucleotide. 
Membrane 92.1 43.8 And the nucleotide homology (91.2%) between IBV 
Nucleocapsid 94.0 33.3 Beijing isolate and IBV LX4 (Accession: AY223860) will 


a) AIBV strain Beaudette (Accession: NC_001451.1); b) SARS decline after the complete genome of IBV LX4 finished. 
eotonavinis dn eressiom NC AOr ee LOn: For there is a gap of more than 6 kb at the 3’end of the 
(iii) Phylogenetic analyses. To determine the rela- | genome, which contains three proteins: S, M and E. These 
tionship among IBV Beijing isolate, SARS-CoV and other sequences have low conservation. The number of amino 
coronaviruses, we compared the overall genome sequence acid of RNA polymerase and four structural proteins was 


HCoV-229E HCoV-299E 
PEDV 
TGEV 
AIBV 
BCoV 
stl MHV SARS-CoV 
(a) 
PEDV. 
MHV RtCoV eee 


HCoV-229E 


HCoV-229E 


Fig. 3. Phylogenetic analyses of coronavirus complete genome and three proteins. Phylogenetic trees were generated by clustw 1.81 and were drawn 
with treeview. (a) Complete genome sequence. HCoV-229E (human coronavirus 229E), NC_002645.1; PEDV (porcine epidemic diarrhea virus), 
NC_003436.1; TGEV (transmissible gastroenteritis virus), NC_002306.2; BcoV (bovine coronavirus), NC_003045.1; MHV (murine hepatitis virus), 
NC_001846.1; SARS-CoV (SARS coronavirus), NC_004718.1; (b) S protein. HCoV-229 E(human coronavirus 229E), NP_073551.1; PEDV (porcine 
epidemic diarrhea virus), NP_598310.1; TGEV (transmissible gastroenteritis virus), NP_058424.1; BcoV (bovine coronavirus), NP_150077.1; MHV 
(murine hepatitis virus), NP_045300.1; SARS-CoV (SARS coronavirus), NP_808851.1; FCoV (feline coronavirus), CAA56778; RtCoV (rat sialoda- 
cryoadenitis coronavirus), AAF97738; CcoV (canine coronavirus), CAA543345; (c) M protein. HCoV-229E (human coronavirus 229E), NP_073555.1; 
PEDV (porcine epidemic diarrhea virus), NP_598313.1; TGEV (transmissible gastroenteritis virus), NP_058427.1; BcoV (bovine coronavirus), 
NP_150082.1; MHV (murine hepatitis virus), NP_045301.1; SARS-CoV (SARS coronavirus), NP_808855.1; FcoV (feline coronavirus), CAA39851; 
RtCoV (rat sialodacryoadenitis coronavirus), AAF97742; CcoV (canine coronavirus), AA033711; (d) N protein. HCoV-229E (human coronavirus 
229E), NP_073556.1; PEDV (porcine epidemic diarrhea virus), NP_598314.1; TGEV (transmissible gastroenteritis virus), NP_058428.1; BcoV (bovine 
coronavirus), NP_150083.1; MHV (murine hepatitis virus), NP_045302.1; SARS-CoV (SARS coronavirus), NP_808858.1; FcoV (feline coronavirus), 
CAA74230; RtCoV (rat sialodacryoadenitis coronavirus), AAF97743; CcoV (canine coronavirus), A44056. 
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also observed. The orfla of IBV Beijing isolate was 36 nt 
longer than that of IBV strain Beaudette (Accession: NC_ 
001451.1), which resulted in 12 amino acids more for 
RNA polymerase than that of IBV strain Beaudette. 
Among the four structural proteins, S protein, M protein, 
and N protein of IBV Beijing isolate were longer than 
those of IBV strain Beaudette (7a.a, la.a, la.a, respec- 
tively). The genome sequence of IBV Beijing isolate was 
very different from that of IBV strain Beaudette, and it 
would be very important for the studies of diagnosis and 
diseases control on IBV in China. 

Moreover, the full-length genome of IBV Beijing 
isolate showed 50.8% homology to that of SARS 
coronavirus (Accession: NC_004718.1). The identities of 
five proteins of IBV Beijing isolate were fairly low at the 
amino acid level, compared to those of SARS coronavirus, 
respectively. Phylogenetic analyses also showed that 
Beijing isolate was not resembled closely to SARS coro- 
navirus, which was consistent with other reports!'' '*), 
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